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I. INTRODUCTION 



Recent experiments on liquid He 4 in vycor and other porous media have revealed 
some interesting features of disordered boson systems. For a boson system at low temper- 
ature, strong disorder may destroy superfluidity and transforms it onto a Bose-glass phase. 
While many theoretical efforts have been devoted to this disorder-tuned quantum critical 
phenomenon — ||, some important issues are yet to be explored and to be understood. 
Like it in the fermion system, study of the nature of low energy excitations is extremely 
important to the understanding of the macroscopic properties. In this paper we study the 
properties of low energy excitations in weakly disordered boson systems and their effects on 
superfluidity. This is of interest both to compare to the pure case, and to study the precursor 
to the destruction of superfluidity. Although the experiment on superconductor-insulator 
transition in dirty films |T(| is an extremely interesting possible realization of disordered 
boson systems ||, the present paper will be restricted to the study of dirty bosons with 
short-ranged interaction, which is suited to model systems like He 4 in random media. 

As usual in theoretical physics, one may wish to proceed the study by starting with a 
perturbation theory about some known limiting cases. However, it is clear that the non- 
interacting system is not a suitable starting point, since disorder will cause a (artificial) 
Bose-Einstein condensation of bosons onto the lowest energy, hence localized state. This 
condensate is unstable with respect to arbitrary weak but finite repulsions between bosons. 
Therefore it is essential to include the interaction between bosons to proceed the study 
of the effect of disorder. Since the condensate is no longer uniform in the presence of 
disorder, the self-consistency required for the normal Hartree theory imposes a set of non- 
linear equations on the condensate which are not solvable in general. In this paper, we adopt 
an alternative approach. We start by considering a model of lattice hard-core bosons with 
random potential, introduced by Ma, Halperin and Lee 0: 

n = -t £ bfbj + h.c + - fibfti , (l) 

<i,j> j 
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where bj, bj are boson creation and annihilation operators at lattice site j, with the hard- 
core constraint bj~bj = 0, 1. < i, j > indicates nearest neighbor and Wj is the random on-site 
potential obeys certain (independent) distribution. This model is equivalent to a quantum 
spin 1/2 XY- model with transverse random fields [0: 

n = - J £ (Sfs* + sfs]) - £ hfl , (2) 

<i,j> j 

with 

Sf <— > bj , 57 <— > bj , 

J < — > 2t , hj < — > fi — Wj . 

Superfluidity in the boson model (|l]) corresponds to a spin long range order in the XY plane. 
Thinking of the boson problem in terms of the spin language and vice versa is proven to be 
a fruitful way in understanding the physics of these problems @,|]||. For the most part of 
the paper, we shall devote ourselves to investigate the ground state properties. 

Above mapping is an exact mathematical transformation. Now we apply approximations 
in term of the spin language. We start by letting the magnitude of the spins to have an 
arbitrary (integer or half integer) value S, roughly corresponding to relaxing the hard core 
condition and letting the site occupation number of bosons bljbj = 0,1,..., 2S. When 
5^00, @ describes a system of classical spins and the usual mean field theory becomes 
exact. We study this Hamiltonian through a 1/5* expansion. When finite S is considered, 
one introduces quantum fluctuations as well as stronger on-site repulsions between bosons. 
To the first order of 1/S, this expansion describes the Gaussian fluctuations and corresponds 
to a generalized spin wave theory, where the magnons correspond to the phonon excitations 
in the boson system. Since in the classical case the long range order persists up to the 
infinitely strong disorder limit (provided that the probability distribution for the random 
field hj, P(hj), is finite at its mean value) 0, the expansion is assumed to be inside the 
superfluid phase. In the end, we will discuss the possibility of probing the transition to the 
disordered phase within this approach. 



3 



Our present study is confined within the first order 1/S expansion, i.e., the spin wave 
theory. In this scheme, the system is described by a quadratic Hamiltonian of (the Holstein- 
Primakoff ) bosons with zero chemical potential. For the pure system, which will be discussed 
in detail in the following section, this theory gives the same results as those in the Bogolui- 
bov theory for weakly interacting bosons. We study the disordered system in section III. In 
this case, the random transverse fields {hj} result in both diagonal and off-diagonal disorder 
for the Holstein-Primakoff bosons. For weak randomness, where perturbation theory is ap- 
plicable, modifications to the low energy spectrum due to disorder are obtained analytically. 
We find that in this regime the superfluidity is rather robust against (weak) disorder, as one 
expects on the physical ground. We conclude our discussions in section IV. 



II. PURE SYSTEM 

In this section we present the spin wave theory results for the pure system. Consider the 
quantum XY-model (|2|) without fields 

n° = - j £ (s*s* + s?s*) , (3) 

<i,j> 

where we have made a spin rotation for convenience. Assuming the ground state to be 
ferromagnetic with broken symmetry in the ^-direction, the classical ground state is then 
given by < Sf >— S. Fluctuations about the classical ground state can be studied through 
the usual Holstein-Primakoff transformation [IT], given by 



S+ = V2S(1 - a] aj /2S) 1/2 a 3 , 
ST = y/2Sa){l - a)a 3 /2S) l/2 , 

S] = S-a}a 3 , (4) 

where a], a 3 - are bosonic creation and annihilation operators. Assuming periodic boundary 
condition on a d- dimensional hypercubic lattice, in terms of the Fourier transformed variables 



£e- fe %} , a fc = ^=E^«, , (5) 



we have, with rescaling J — > J/S* 2 , 

H° = - Z -NJ + ^ E 4a* - J£ E 7k(4a* + °fc a -* + H - c -) + ( o§/2 ) , ( 6 ) 

k 1 

where 



7* = (7) 

z s 

S above is the directional unit vector of the lattice, and z is the lattice coordination number. 
Ignoring the terms of higher power in l/S, the usual Bogoluibov transformation gives 

n° = E° + ^ k bib k , (8) 

k 

with 

^ = ^E(0^-i)-f^ > ( 9 ) 

zJ 



^ = T v /1 -^- (10) 

Here b k , b k are the new creation and annihilation operators respectively, with 

b k = u k a k - v k at k , bl = -v k a- k + u k at , (11) 

and 

"2 =5(7^) • ^(tt^H • 2 — Tr^T- (12 » 

The low energy excitations are magnetic phonons (magnons) with speed of sound c° = 

To have some physical feelings about the nature of the low energy excitations, we express 
the magnon operators in terms of the spin operators (to the leading order of l/S): 

b * = ^ " Tfc)- 1/4 SZ " (1 " 7*) 1/4 S£) , (13) 
bt = -L (-<(! - 7jfc )-V^ - (1 - 7fe )V4^) > (14) 



where Sk is the Fourier transform of Sj. S% and S\ generate out-plane oscillations and 
in-plane rotations respectively. Above decomposition of the magnon operators shows that 
while high energy excitations are composed of both kind spin motions, the low energy, long 
wavelength excitations are dominantly in-plane rotations. The zero mode {uj® =0 = 0) is a 
pure in-plane uniform rotation. As we will see later, such a zero mode also exists in the 
disordered system. 

To calculate the helicity modulus, which is proportional to the superffuid density |fP2| , 



we impose a phase twist on the order parameter. Let us maintain the periodic boundary 
condition in the first d — 1 directions but impose the anti-periodic boundary condition in 
the d-th direction for the in-plane component of the spins. Since the spin wave excitations 
are well defined only with the correct (classical) ground state, we need to determine first the 
classical ground state under such a boundary condition. It is clear that the spin configura- 
tion which gives the lowest energy under such circumstance is the state in which spins are 
rotated gradually from angle to ir in the d-th direction (See the next section for discussions 
for general situations.). Under such a rotation, spins in each d — 1 dimensional plane per- 
pendicular to the axis of the phase twist are rotated by the same amount, and the angular 
difference of the rotations between two successive planes is vr/L, where L is the linear size 
of the lattice. To apply the Holstein-Primakoff transformation (f|), we perform therefore a 
spin rotation around the y-axis (in the spin space), and the angle of rotation on the spin at 
site j = (ji, ....,jd-i,jd) is 

= *k . (is. 

In terms of the rotated spins, @ with anti-periodic boundary condition becomes 

n° = -JJ2 [cosih - h^s* + sfsf) - « n fa - h)(s;s; - s*s;)} . (ie) 

<i,3> 

The spin variables in ( |T6D satisfy the periodic boundary condition, and they are ordered in 
the ^-direction in the ground state. 

Applying the Holstein-Primakoff transformation (^), and through a procedure similar to 
that in the case of the periodic boundary condition, we have, to the order 1/S, 



H'^E'V + ^'Xb, , (17) 

k 



with 



J 



25 „ 



u'l = - 2 + 2a)[(z - 2 + 2a) - (z - 2)^ - 2a 7 I] , (19) 

« = «»(") , 7l = 5 £ > ^ = ^ £ e<w± ' 

<5|| 5x 

and ||, _L are defined with respect to the axis of phase twist. 

Compare above results with those obtained with the periodic boundary condition, we 
see that, by twisting the phase of the order parameter, one lifts the ground state energy of 
the system but lowers the excitation energy (with respect to the new ground state energy 
level), i.e., E'° > E°, u'® < u k . Since to the first order of 1/S the second term in (|i6|) 
vanishes, we can think of the effect of changing the boundary condition as a reduction of 
the spin coupling in the d-th direction, i.e., J — > J in the first d — 1 directions and J — > a J 
in the d-th. direction. Because of this reduction, one loses ground state energy but makes 
excitations easier. Moreover, since this coupling reduction is anisotropic, the system picks 
up an easy direction of excitation along the d-axis. 



E = 7^ £( V (* - 2 + 2a)[(z -2 + 2a)-(z- 2)^ - 2a 7 



{z-2 + 2a)) - Z ~y L NJ - aNJ , (18) 



The helicity modulus 7 ||T2[ can now be easily calculated from the expressions obtained 
above: 

7 (T) = lim 2(-) 2 — = lim 2(-) 2 - ( AE° + [ d d k ) , (20) 

where AF refers to the change of free energy as a result of changing boundary condition 
and AE° = E'° - E°, Auj° k = lu'I - uj° k oc (f ) 2 k. This gives 

7 (T) = 7 (0) - aT d+1 , (21) 
where a is a positive constant. This is in agreement with the known result for weakly 



interacting boson systems |13 



Before we go further, we would like to comment on the validity of the spin wave approach 
to our problem. In one dimension (ID), the quantum correction to the in-plane magnetiza- 
tion, 5m = < a ] a j >> is divergent [plf] , reflecting the lack of true long range order in 
the quantum XY model jllj. For two and higher dimensions, it has been shown rigorously 
that finite magnetization exists in the XY plane ||16|| , and indeed the quantum correction to 
the classical magnetization in the spin wave theory becomes finite also. The fact that 5m 
decreases with increasing dimensionality indicates that spin wave theory becomes a better 
approximation in higher dimensions. Notice that despite the divergence in 5m, and thus 
the absence of the condensate, the superfluidity that exists (at zero temperature) in one 
dimension manifests itself through the linear behavior in low energy excitations. 

Finally, we give the spin wave solution for the XY model in a uniform transverse field 
where hj = h^. This can be obtained in a similar manner to the procedure for the case of 
ho = shown above, with a rescaling J — > J/S 2 , h — > h /S. Comparison of the results 
obtained here with those calculated in the presence of the random fields to be described in 
the next section will help us isolate the effect of disorder. For large field, i.e., for ho > zJ, 
the classical spins are completely aligned with the transverse field and therefore there is 
no in-plane magnetization. For weaker field, the classical ground state (with the periodic 
boundary condition) is given by < S? >= ScosOq, < Sj >= SsinO^, where the transverse 
field h tilts the (classical) spins with an angle 8q such that sin8 = H = h Q /zJ. The 
ground state energy 

The magnon spectrum is modified by ho as 

zJ 



ul{h a ) = -^l- lk ){l- lk H%) . (23) 

One can compute the out-plane susceptibility x± (compressibility in the boson language) by 
taking second derivative of E°: 

f)2pO/Tvr C2 



The S 2 term is simply the classical result, and the rest are corrections due to quantum 
fluctuations. Since the classical contribution to 7(0) is just J, we calculate the speed of 
sound 

< ~ x± ~ & , (25) 

in agreement with the result obtained earlier. We see that from pure classical considerations 

one can obtain c s which characterizes the quantum nature of the low energy excitations. 

Not surprisingly the speed of sound is reduced by h , with c° s (h ) = c° s ^Jl — Hq. However, 

it is less obvious that the (relative) quantum correction to the in-plane magnetization 

Sm(h ) 1 _ t 1 f d d k ( l-fr(l + flg) \ 
— E < a l a k >= 7^ / 77^ 1 = ~ 1 > ( 26 ) 



Scos6 NS^ * K 2SJ (2n)« \ ^(l_ 7fc) (l- 7fc ^ 

is also reduced, i.e., < ^51. This suggests that the reduction of < S? > by a uniform 

transverse field is mainly a classical effect, resulted from the tilting of spins off the ZX (XY 
in the original coordinate system) plane by ho. In fact, as H — > 1, so that cos9 — > 0, 
5m/ ScosOq — > 0, which helps to explain why the critical value of H is just the classical 
value. We will compare this with the effect of random fields as we proceed. 

In summary, the quantum XY model provides us a useful representation of the boson 
problem in studying the physics of superfluidity. In the large S theory, the expansion 
parameter 1/S contains both the effect of the repulsive interactions between bosons and 
the effect of quantum zero point fluctuations. We have shown that for the pure system the 
quadratic fluctuations represented in the first order 1/S (spin wave) theory gives the same 
physics as that described in the Bogoluibov theory for weakly interacting bosons. Now we 
apply this method to study the disordered case. 

III. DISORDERED SYSTEM 

Consider the disordered system (with periodic boundary condition) given by 

w = -i E ( s ' s i + s i s !) - £ tH* > (27) 

15 <i,j> j D 
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with 

Tj = , ~h~h~j = 5ijh 2 . (28) 

In (PTD, we have made the rescaling J — ► J/S 2 , hj — > hj/S. The overline indicates average 
over the random transverse fields. The choice of ft^ = corresponds to a 'particle- hole sym- 
metry' in the boson system. For hard-core bosons, it corresponds to an average occupancy 
of half particle per site. Assume that spins in the ground state of the pure (classical) system 
are ordered in the z-direction, the transverse random fields along the |/-axis then tilt each 
(classical) spin according to the equation of motion (with §f assical = S(0, sin9j, cosOj) ): 

sinOjJ ^ cosOji = hjcosOj , (29) 

<?> 

where < j' > indicates nearest neighbors of the site j. Now we perform a local rotation of 
the (quantum) spin Sj over an angle 9j about the x-axis so that the classical ground state 
is uniformly ordered in the z-axis in terms of the rotated spins. After the rotation fl27|) 
becomes 

H = -lL E ( S i S l + cosOiCosOjS'S* + sinOiSinejSfS] 

S <i,3> 

-sine.cosejSfS] - cosOisinO^SfS]) - ]T MsinOjS* + cosOjS]) . (30) 

3 S 

Applying the Holstein-Primakoff transformation (f|) now, together with (|29|), one has 

7i = — J cosOiCosOj — ^2 hjSinOj — ^— ^ [(1 — sin9isin9j)ai(ij + 

<i,j> j <M> 



1 h ■ 1 

(1 + sinOisinO^alcLj + H.c] + - ~f a ) a i + ( 0375) • ( 31 ) 

j •* 

Rewrite ([!!]) in terms of Fourier transformed variables, we have 

n = n° + n 1 , (32) 

with H° given by (^) and 

U l = E 10 + Y1 u k - k ia\a k i + Y^{ v -k-k'a k a k > - v k ^ k >a\a k i + H.c) , (33) 



h h' h h' 
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where 



E 10 = J ^2 (1 — cosQiCosOj) — hj-sinOj , 

<i,j> j 



1 Vf— j 



NS ^- jK sinO j 

3 J 



zJ)e 



'/h.-.L > 



Vk,k' 



J 



2NS 



£ sin9 i sin9 j e i( - kXi+ ^ . 



(34) 



(35) 



(36) 



<i,j> 



In the weak disorder limit, H 1 may be considered as a perturbation to TC°. 

In stead of carrying out a Bogoluibov type transformation, which becomes non-local and 
depends on the configuration of the random fields in the disordered case, we compute poles 



of the propagator < T[ak(t)al(0)] >(<...> denotes the ground state expectation value 
and T indicates time ordering.). These poles are directly related to the excitation spectrum 
of the system. To this end, we consider 



F(k, k';t) = -i 

which obeys the equation of motion 

,.d 



1 <T[a k (t)at,(0)}> ' 



< T[a+^)a+(0)] > 



(37) 



i— - T k )F{k, k'; t) = 5{t)5 kk ,E l + £ U kr F(k" , k'; t) 
at k „ 



(38) 



where 



T, = 



zJ 



^ 1 7fe 


Ik \ 




( A 


2 


2 






Ik 

\ 2 


-a-?)y 







(39) 



and 



U k h' 



Uk-k' — Vk-k' — V~k',k V k - k < + V-k',k 

y —(Vk,-k' + V-k',k) —{Uk-k' — V k - k > — V- k ',k) j 

Applying the Fourier transformation in time, with 



(40) 



F(h,k';u) 



dtF(k, k'\ t)e iuJt 



(41) 
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we have 



or 



(u - T k )F(k, k'\ u) = 5 kk ,E l + U k , k »F(k n , k'; u) 

k" 



F(k, k>- u) = {u- T k )-H kk ,E l + T,)- 1 ]T U kr F(k" , k'\ 



k" 



For the pure system, U kjk i = 0, we have 



(42) 



(43) 



F u (fc, k ; u) = F u (fc, u))5 h ,k>Ex = G u (k, u) 



(44) 



with 



F°(k,u) 



T k )~ l 



D°(k,u) 



zJ 7k 
S 2 



and 



zJ 7fc 
S 2 

S v 1 2 J 



(45) 



D°(k,u) (u-u% + iri){u + ul-iri) ' 



This gives 



G°(k,u) 



(46) 



(47) 



y UkVk \ w-ul+ir! U+cJ^in J J 

where u k and v k are given by fll2|) . Thus, from the poles in G°(k,cj), we have recovered the 
excitation spectrum obtained earlier through the Bogoluibov transformation. 

For disordered system, one can iterate (f43|) to find a perturbative solution for F(k, k'\ iS). 
U kjk r contains randomness that needs to be averaged out. It depends on the random fields 
{hj} both explicitly and implicitly in {dj} through (|29|). For weak randomness, where 
h « zJ, one can solve ( p9|) order by order (see the appendix). One has, in the regime of 
weak disorder, 



G(k, u) = F(k, k- u) = (F°-\k, lu) - E*{k, uo)Y 1 E 1 



(48) 
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where the self energy H*(k,u) is given by 



Z*(k,uj) = U k , k + J2Uk,k'F (k>,Lu)U k ,, k . (49) 
k> 

The detailed algebraic expression of T,*(k,u) is given in the appendix. This is consistent 
with a perturbation expansion for the self energy to the order h 4 (see the appendix). Notice 
that in real space U k ^ is not just a random on-site potential. It contains both (correlated) 
diagonal and off-diagonal disorder. 

The poles of G(k,u) are given by ( |A15| ) when the distribution of the random fields is 
Gaussian, which can be computed analytically for ID system. As discussed earlier, while 
the divergence in 5m indicates the instability of the classical ground state and thus the 
absence of true long range order in the quantum system, the linear behavior in uj® at small 
k indicates the existence of the superfluidity at zero temperature. It is therefore meaningful 
to investigate deviations from u) k due to disorder in ID systems within the present approach. 
To the order h 4 , the poles of G(k,u) are given by (with z = 2) 

u = u h = ul[l- l -{^YA{k)\ , (50) 

with 

Re{A} = ^ + 5cosk - 3cos 3 k + cos 4 k , (51) 

T r .-i (2 — cosk) 2 . . ... , . .„ 9 ,.l — cosk . 

Im{A} = - [sink (1 + cosk) + (1 + cos 2 k)— — . (52) 

2 A\smk\ 

We see that the linear mode persists at low energy with a reduced speed of sound c s = 
c°(l — y (tj) 4 )- It has been shown |4]] that for S = 1/2 any amount of disorder will change the 
power law behavior of the spin-spin correlation function to an exponential one, corresponding 
to the instability of superfluidity in the ID hard-core boson system. For systems with soft- 
core bosons, which are roughly described by (0) with S > 1/2, renormalization group study 
T7J shows that superfluidity may persist in the presence of (weak) disorder, if the value of 



an exponent rj of the correlation function in the corresponding pure system is less than a 
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critical value r\ c = 1/3. The fact that in our calculations the linear mode persists in the 
weak disorder limit suggests that present 1/S theory describes systems with S > S c such 
that the exponent r\ < r\ c . The imaginary part of the pole diverges at the zone boundary 
(k = ±7r), indicating vanishing life time of these (Bloch) modes. By analyzing the scattering 
rate of the Bloch phonons using Fermi's golden rule, this divergence can be interpreted as 
the divergence of the density of states in one dimension. In 2D, the imaginary part has a 
(van Hove) singularity at the band center, and it behaves regularly in higher dimensions. 
For an arbitrary dimension d, the correction to the speed of sound is given by 



* = * H^H • 



(53) 



with 



m = ^l +V l^L ^^l . (54) 
K ' 2z J {2n) d l- 7fc V ; 

The imaginary part of D(k,uj) = de£(F 0_1 — £*) vanishes for small k and uo as uo d+2 (see 
the appendix), implying a decay rate of phonons r _1 oc u d+1 , and the mean free path 
I = c s r oc k~( d+1 \ Thus for small fc, kl » 1. This implies that momentum is still 
approximately conserved. We see that both the reduction of the sound speed and the decay 
of phonons become weaker in higher dimensions. Notice that this is not a time decay rate 
for phonons, which can only be nonzero through multi-phonon scattering, but is just the 
scattering rate for a Bloch phonon. Since the scattering is elastic, energy conservation 
implies r _1 oc a/* -1 , the extra suppression factor of uo 2 is presumably a reflection that there 
always exists a zero frequency mode which corresponds to uniform rotation in the ordering 
plane. Since for small u>, ujt » 1, the Bloch modes remain robust in the presence of weak 
disorder. The fact that the zero mode (oUk=o = 0) has zero imaginary part shows that it 
remains to be an exact eigenmode of the Hamiltonian in the disordered system. 

Next we study the effect of disorder on the magnitude of the order parameter, which 
corresponds to the square root of the condensate density in the boson language. It is given 
by 
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m = < S* > = j-j:'^2cos6j(S— < djdj >) , (55) 

3 

which can simply be understood as following: the first term gives the reduction in m due 
to the tilt of the classical spins from the ordering plane. The second term shows that each 
Holstein-Primakoff boson lowers the spin along z'-axis by one, and hence only by cosdj along 
the z-axis. Thus to the order h 2 

m = —Y, cos6j(S -j7^2< 4°* >o) 

3 k 

where < a\ak >o is the expectation value with respect to H° (since 5 < a\a k > vanishes to 
order h 2 , see ( [MT]) .). 

Again let m = m dassica i — 5m, where m classica i is the classical (S — ► oo) condensate, 
and 5m is the reduction due to quantum fluctuations. In the pure case where hj = h and 
9j = 6q, we have seen that besides the 'classical' factor cos8q, 5m is also reduced by Hq, 
resulting in a net reduction of 5m/ mdassicai by a uniform field. This is in contrast to that 
in the random system considered here, where 5m/m c i assica i is independent of the disorder to 
order h 2 . While 5m is reduced by disorder to order h 2 , indicating that disorder and quantum 
fluctuations have opposing effects on the condensate density, it will not be the case to higher 
order in h, and that 5m will also be enhanced by disorder. To order /i 4 , we have (see the 
appendix) 

5m 1 r d d k ( 1 - ^ \ 1 h 4 r d d k 1 

m clasmcal ~2SJ (27p WT^Tk J + 4S ( i7 } J (27c) d V^Yk 

l^(T(1r\ T(1r\-L l T(V\\ (l-f) 2 (l + / 2 (fc)) f ~ \ 

{^{H k )-h{k) + -h{k)) \- lk + 1 ~ V 1 -7kh{k) J> , (57) 

where I\(k), ^(k) and ^(k) are given in the appendix. The fact that the second term in 
the right hand side of (|57P is positive suggests that disorder enhances quantum fluctuations. 

One can calculate the helicity modulus in a similar way as one does for the pure model. 
In the disordered case, the effective Hamiltonian for phonons become considerably more 
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complicated for a system with anti-periodic boundary condition. For classical spins, the 
equations for the ground state are given by 



<j'> 



sinOj ^2 cos(pjicos6ji = -j-cos(J)jCos8j , 



sincf)j 2J coscj)j'Cos9j> = coscpj 2J sin<fij>cos9ji , (58) 
<f> <j'> 

where the classical spin is defined by Sf assical = S(sin(f)jCos9j, sin9j, coscj)jCos9j). Obviously, 
we expect the phase twist to be large where the tilt of the spins away from the ordering 
plane is large and vice versa. For a general random field distribution, ( |58"D may be solved 
numerically. The helicity modulus 

7 (T = 0) = J(l - a (±)* + Q{h')) , (59) 

with a > 0. In the special case where {hj} is given by a bimodal distribution, P{hj) = 
h(S(hj — ho) + 8(hj + ho)), the solution for in (^) is the same as that for a uniform 
field ho, and a = 1. Since there is no shift in the phonon spectrum to the order h 2 , there 
is no quantum correction to 7 to that order in the disorder. One might ask if the (low) 
temperature dependence of 7 is affected by the disorder. However, since the low energy 
excitation remains to be phonon like, the density of state N(u) oc uj^ 1 , assuming that the 
single mode approximation remains intact. This implies (see (PH)) that the T d+l behavior 
of the temperature dependence of the helicity modulus should be unaltered by the presence 
of weak disorder. 

To summarize our results, we find that to the lowest non- vanishing order in h/zJ, and 
relative to the pure system, the speed of sound is reduced (5c s oc (jj) 4 ), the condensate 
density is reduced classically, but unaltered quantum mechanically, while the superfluid 
density is reduced classically, but remains unchanged quantum mechanically. To higher 
order of h, disorder tends to enhance quantum fluctuations. 

So far we considered only a special case where hj = h = 0. In general, h is nonzero 
and instead of (|A1|) , we have 



16 



% = ^E^^-4^ + ^ 2 ) > (6°) 



with 



Compare with (AT), we see that the primary role of a non- vanishing /i is to introduce 



correlations among 0j's. Calculation with finite /io is more complicated. To see physics, 
however, one can consider the case of correlated randomness. Thus, we consider (p8|) with 



hj = , Zii/ij = f(\xi - xj\) . (61) 



where f(x) is some arbitrary function. To the lowest order, poles of the Green function ([481) 
are given by 

^ 2 = (^) 2 (l-7 fc )(l-7 fc (^) 2 ) , (62) 



where h 2 is the average of a pair of random field at nearest neighbor, hjhj + s = h 2 . Thus, 
in the weak randomness limit, correlated random fields reduce the excitation energy, and 
hence the speed of sound in the same manner as a uniform field (see (|23|)), and do not induce 
finite scattering rate for the Bloch phonons at the lowest order. (Note that in ( |62| ) only the 
nearest neighbor correlation in hj matters.) 

IV. DISCUSSIONS 

As remarked previously, our approximation consists of a double expansion in 1/S and 
the strength of disorder h about a saddle point solution of ( |27"D which becomes exact in the 



classical (S = oo) limit. Now we re-examine this approximation scheme in terms of the 
original boson Hamiltonian (|TJ). In terms of the hard-core boson operators, the classical spin 



ground state described by (p9|) corresponds to a Gutzwiller-type trial wavefunction |18 



>= Y[(simpj + cosLpjb+)\0 > , (63) 
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with (fj = it/A — 6j/2. Minimizing the energy E =< ^/g\T-C\^/g > with respect to {9j}, one 
recovers (|29|). In this state, the order parameter 

6 = 4<*gIE^I*g>= 2^E«>^ > ( 64 ) 

which is equivalent to the in-plane magnetization of the classical spins (with S — 1/2) con- 
sidered previously. The quantum fluctuations described by the spin-wave theory correspond 
to the Gaussian fluctuations about this Gutzwiller state. Within this scheme, our pertur- 
bative calculation shows that the low energy Bloch modes are rather robust, and disorder 
has little effect on the long wavelength quantum mechanical behavior of the system. We 
cannot, however, rule out the possibility of a stronger effect of disorder (e.g. a h 2 correc- 
tion to the speed of sound) as one goes to higher orders of 1/S. At higher orders of l/S, 
phonon-phonon scatterings take place and cause decay of these quasi-particles. Presumably 
such decay processes are stronger than those of the Bloch phonons in the corresponding pure 
system due to disorder enhanced quantum fluctuations. How does disorder enhance these 
fluctuation effects is an important question that yet to be investigated. 

Our perturbative study shows that superfluidity is rather insensitive to the presence of 
weak disorder. This result is fully in agreement with what one would expect intuitively: 
since superfluidity is due to quantum coherence on the macroscopic scale, superfluid should 
be rigid against weak impurity scatterings. A similar statement for fermion systems is 
expressed in the content of the Anderson Theorem for the BCS superconductors with weak 



homogeneous non- magnetic impurities fO^ . A natural question to be asked is then how 
would the system evolves with increasing randomness. At zero temperature, as disorder 
becomes stronger, or as effects of interaction and quantum zero point fluctuations become 
more important, or both, one may reach a point where a transition from the superfluid 
phase to a disordered (Bose glass) phase takes place. Since in the Bose glass phase the low 
energy excitations are single particle like, one may expect that the speed of sound, which 
characterizes the low energy excitations of the superfluid phase, vanishes at that point. Thus 



the transition point can be located naively by setting c s in (|53| ) to zero, which gives 



\^Y-m = i . (65) 

Of course, long before this point is reached our (weak disorder) approximation loses its legit- 
imacy. However, it is still a non-trivial and interesting question that whether the Gaussian 
fluctuations contained in the first order 1/S expansion are capable to describe such a phase 
transition. Since this approach is based upon an expansion about the ordered phase, phase 
transition is expected to be signaled by instabilities of the expansion, such as appearance of 
negative energy modes or strong divergence of quantum zero point fluctuation corrections. 
The present weak disorder calculation cannot answer such questions concerning strong dis- 
order, but this approach does provide a scheme for further investigations. Since the effective 
Hamiltonian is quadratic to the first order of 1/S, exact numerical diagonalizations are pos- 
sible for finite systems upto sizes which are unreachable otherwise. Thus it provides us a way 
of study the low energy excitation spectrum in the strong disordered system. This will be 
discussed in a forth coming work [pOfl . Of course, such a study can only address these ques- 



tions within the Gaussian scheme. Other than simply having a vanishing speed of sound, 
there are different possibilities for the phonon mode to evolve into the single particle con- 
tinuum as the system is tuned into the disordered phase. For instance, c s may remain finite 
while phonon decays strongly with increasing disorder so that its spectral weight vanishes 
at the transition point. How precisely the phonon mode evolves with increasing random- 
ness is still an open question which can only be answered by going beyond the Gaussian 
approximation. 

In conclusion, we have found that the spin representation of the boson problem is an ef- 
fective approach for investigating the effect of disorder. Within the Gaussian approximation, 
our perturbative study shows that the superfluidity remains robust in the presence of weak 
disorder, while random fields have scatterings with Bloch phonons with a rate proportional 
to u d+1 and give a weak reduction of the speed of sound. Our calculations also suggest that 
stronger disorder tends to enhance quantum fluctuations which may eventually cause the 
destruction of superfluidity 0. 
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APPENDIX A 



In this appendix we evaluate the self energy E*(fc,u;) to order h 4 . Assuming h is small, 
one can solve (18) perturbatively. To the lowest order of hj, 9j is simply given by Hj = hj/ zJ . 
To the next order, one has 



' <j r > 



(Al) 



Here we see that Oj starts to couple to the field at the neighbor sites as one goes to higher 
order of hj. We calculate the quantities 



^=^(# 4 -(# 2 ) 2 ) , 



(A2) 



TC=(^) 2 ^(l-7 fc ) 2 (^ 4 -(^ 2 ) 2 ) , 



(A3) 



(v k - k > + v_ k ^ k ){v k ,_ k + v_ Kk >) 



:f) 2 i ^(i+7, + ,)(/? 5 ) 2 , 



where H n is the n-th moment of Hj. With these results, the self energy 



is evaluated as 



Z ^(w- {Wf) + (g) 2 ^(^ - m 2 ) x> + f a - f )) 

(l- 7fc _ fc /) 2 zJ 2 2 Z J 1 + 7 fc _y 

D°(k',u) l + [ 2S } zN [ } % S [ lk,) D°(k>,uj) ' 



(A4) 



(A5) 



(A6) 
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= — - (^) 2 ) + (g) 2 ^ - (^) 2 ) H(" - ^(i - f )) 

(l-7 fc - fc 2 zJ 2 2 Z J 1 + 7fc _ fc / 

£>°(A;',w) J l 25 J ^iV l J ^ S l D°(k',uj) ' 



(A7) 



(-) - ) ) L IY D0(fc>) 
-(g) 2 ^(^) 2 E^(i-7 fc 1 + 7fe - fc ' 



if s 



with D°(k,uj) given by fl46[) . Since 

/ 



D°(k',co) 



_( zJTfe _ sp* \ 
\ S 2 ^12^ 



the poles in G(k,u) is given by the zero's in 



V 



+ E^) W _f(i-^)_E^ 



(A8) 



(A9) 



j D(A;,cj) = det(F 



o-i 



E*) = u/ 



1_ 5<if ) 2 (H*-m 2 )i2(k,u>) 



-(f ) 2 (1 - 7.) - i(f ) 2 (^- (^) 2 )(1 - |)(1 + (f ) 2 / 2 (A^)) 
+ 1 {Z 4 {W)2){1 ~ 7*)Js(M) - ^(f) 4 m 2 (l - 7 fc )/i(^^) , 



(A10) 



with 



r/. \ 1 V- (1 -7fe0(l +7fc~fc0 

Jl(M = ^ ^(^0 



(All) 



(A12) 



r ,, s _ 1 ^ iv (1 - 7fc-fcQ 2 
When the distribution of random fields is Gaussian, 



(A13) 



H* - (F 2 ) 2 = 2(-^-) 4 . 

/6 t/ 



(A14) 



To the order 0(h 4 ), one can substitute uj in Jj (j = 1,2,3) by uj®. Since we are only 
interested in the case where uj > 0, we may ignore the singularities for negative uj. We have 
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UJ = UJ 



02 



1 + 



zJ' 



h{k) + (1 



7fc 1 + / 2 (fc) 
2 J 1-7* 



7 3 (A;) + 



where (rj — > 0~ 



/i(fe)- 1 E (1 ~ Tfc,)(1+Tfe ' fc0 



(A15) 



(A16) 



i _ (i - 7 *-*0 2 



N Tf Ik' - Ik + ir] 



(A17) 



/ f m = — v 7fe/ ^ - 7fc - fc ') 2 

31 J ^ W ^ 7*' — 7* + ^7 



(A18) 



To the order of a; 2 and A; 2 , 



^ 2 = a;° 2 1 



/i 4 3z + 3 1 x . Ej=i «n 2 fc; 



,21,/ 



2^ 



iV^ 1 



(A19) 



fc , ^ 7 fc / - z?7 7 

which implies a correction to the speed of sound due to disorder shown in eqn. fl53|). The 
imaginary part of ( A10| ) is of the order u d+2 , which can be seen through a simple dimensional 
analysis. This gives the life time of phonons oc u~( d+1 \ 

To study the effect of disorder on the quantum corrections to the order parameter 
one needs to compute 

I 

' k 



k 



(A20) 



Using the result for G obtained above, and completing the ui integral by the usual contour 
integration over the semicircle at the lower half complex plane, we obtain 

I p (jfik f 1 — 241 ^ 1/1 ' Adl 



1-1 \ 1 h 4 r d d k 1 

< aJaj> ~2J (2ny[^T=T k 1 ) + 4 { zJ } J (2vr)^ V^Tk 



l-lkh(k) 



(A21) 



1-7* 

One can verify numerically that the second term is positive. From this result, it is easy to 
show that the relative quantum correction to the order parameter is given by (|57Q . 
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